Multipartite entanglement detection in bosons 
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We propose a simple quantum network to detect multipartite entangled states of bosons, and 
show how to implement this network for neutral atoms stored in an optical lattice. We investigate 
the special properties of cluster states, multipartite entangled states and superpositions of distinct 
macroscopic quantum states that can be identified by the network. 
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Quantum entanglement is a very important physical 
resource in quantum information processing tasks, such 
as quantum cryptography teleportation, quantum fre- 
quency standards or quantum-enhanced positioning 
Since the implementation of any of these tasks requires 
precise knowledge on the entangled states being used, the 
development of "measurement tools" for the characteri- 
zation and detection of entanglement in physical systems 
is of great practical importance. 

While the creation of multipartite entanglement has 
been recently achieved in a controlled way in experiments 
with either Mott insulating states of neutral bosonic 
atoms in optical lattices @> 0] or ion traps 0], its un- 
ambiguous detection has so far not been possible @ in 
these setups. The usual experimental methods to detect 
entanglement are based on the violation of Bell type in- 
equalities which are known to be quite inefficient, 
in the sense that they leave many entangled states un- 
detected 0. Alternatively, one can perform a complete 
state tomography of the system Q , but this method re- 
quires the preparation of an exponentially large number 
of copies of the state and it is redundant, since not all 
parameters of the density operator are relevant for the 
entanglement detection. 

In this Letter we propose a simple quantum network to 
detect multipartite entangled states through an entangle- 
ment test more powerful than the Bcll-CHSH inequalities 
for all possible settings ||, albeit less powerful than full 
state tomography. The network is realized by coupling 
two identically prepared ID rows of N previously entan- 
gled qubits via pairwise beam splitters (BS) , as shown in 
Fig. H We also show how to implement this network in 
an optical lattice or array of magnetic microtraps loaded 
with atoms in a Mott insulating state with filling factor 
one ^3 • Each of the atoms has two long lived internal 
states a and b which represent the qubit. The pairwise BS 
can be implemented by decreasing the horizontal barrier 
between the two rows of atoms. 

We start by introducing the set of inequalities used by 
our network for the detection of multipartite entangle- 
ment. The information-theoretic approach to separabil- 
ity of bipartite quantum systems leads to a set of entropic 
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FIG. 1: Network of BS acting on pairs of identical bosons. 
The two rows of N atoms, labelled I and II respectively, are 
identical, and the state of each of the rows is pi23...N- The 
total state of the system is P123...N ® P\23...n- 



inequalities satisfied by all separable bipartite states |8(. 
We extend these inequalities to separable multipartite 
states by considering a state P123...N of N subsystems. If 
P123...N is separable then we can write it as 



P123...N 
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Cep[ ® P2 ® p\ ® ■ ■ ■ <£> p N , 



(1) 



where p^ is a state of subsystem j, and J2 e Ct — 1. The 
purity Tr(p? 23 „) of pu3...n, where n € 1,2, ...,N, is 
smaller or equal than the purity of any of its reduced 
density operators. For example, 
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Tr (Pl23...n-l) < Tr (Pl23...«-2) 

<Ti-(p 2 12 )<TT(pl)<l. (2) 



This set of nonlinear inequalities provides a set of neces- 
sary conditions for separability, i.e. if for any state p any 
of these inequalities is violated then p is entangled. For 
the case where p\23... n is separable and pure we have that 
Tr (pi23...„) = 1 and the inequalities become equalities. 

In order to test Eq. (J2J we need to be able to determine 
the non- linear functional Tr(p 2 ), where p is any of the 
different reduced density operators of pi23... n - The direct 
estimation of this functional has been addressed in 
where the value of Tr(p 2 ) is determined by measuring 
the expectation value, on state p ® p, of the symmetric 
and antisymmetric projectors. 
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Let us first consider the simple scenario of one pair 
of identical bosons, in state pj ® pj, impinging on the 
BS (as depicted in Fig. P). After applying the BS we 
write the purification of the state p'j = UbsPj ® Pj'^bS' 
where Ubs is the unitary time evolution operator of the 
BS, in the triplet/singlet basis. The operator Ubs trans- 
forms a(b)^j -» {a(b) ( /}j - ia(b)% \ Z )V2, where a\ j) 

and 6j are the bosonic destruction operators for par- 
ticles in row Z = J, J/, site number j = 1 • • • N and 
internal state a and 6, respectively. The triplet basis 
elements are (a^ )t /3p )t + affi pffi)\v&c), a,/3 e a, 6, 
with respective pre-factors c^, and the singlet element 

is (a^b^ - a^ )+ 4 j)t )l vac >> with pre-factor c^°, and 
| vac) the vacuum state. The symmetric/antisymmetric 
components of pj (g) pj are spanned by the triplet /singlet 
basis elements. Note that in the triplet basis states, 
both bosons are in the same spatial mode, while in the 
singlet basis state each boson occupies a different spa- 
tial mode. Hence, measuring the probability of find- 
ing two bosons in, cither the same spatial mode or dif- 
ferent spatial modes, after the BS, is tantamount to 
measuring the probabilities P± of projecting the state 
Pj pj on its symmetric or antisymmetric subspaces. 
These probabilities are given by the expectation value 
of the symmetric/antisymmetric projectors S± defined 
by S± = (I ± V)/2, where I is the identity opera- 
tor, V is the swap operator, with Vaj /3jj | vac) = 

Pi = ^Tr [(I ± V) Pj ® Pj ] = \ ± ^Tr {p% (3) 

where " — "(" + ") stands for projection on the 
(anti)symmetric subspace. Thus, after we let pj ® pj go 
through the pairwise BS we can determine the purity of 
Pj by measuring the projections of p'j on the symmetric 
and antisymmetric subspaces. For pairs of polarization 
entangled photons, a scenario similar to this one is cur- 
rently being implemented experimentally [Icj . 

We now extend the above two-boson scenario to the 
general situation (see Fig. and consider two copies of 
a state of N bosons, undergoing pairwise BS. By cor- 
relating the probabilities of projecting the state of each 
pair of identical bosons on the symmetric/antisymmetric 
subspaces, we can estimate the purity of P123...N and of 
any of its reduced density operators. As a more concrete 
example, let us consider the probabilities for N = 3. We 
will label the subsystems 1,2,3, respectively: 
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P± 1±2 ± 3 = ^Tr[JJ(I±^-)Pi23®Pi23] 
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i[l± 1 Tr(^)± 2 Tr(^)± 3 Tr(^) 
± 1>2 Tr(p 2 12 )± lt3 Tr(p 2 13 )± 2t3 Tr(pt 3 ) 



(4) 



±i,2, 3 Tt(Pm)L 
(±0(±i), = 



1,2,3, and V\ 2 3 stand 



where ±i,j' 

for the swap operator acting on subsystem 1,2,3. The 
purities related to P123 are unequivocally determined 
by the eight probabilities P± 1 ± 2 ± 3 - The expression for 
the probabilities in Eq. Q can be straightforwardly ex- 
tended to states of N bosons, where we consider the 
expectation values of the projector J^£ i=1 (][ ±^ Vi)/2, on 
P123...N <8> P123...N- In the N boson case, the 2^—1 un- 
known purities will be determined by the 2 N — 1 inde- 
pendent probabilities. 

The implementation of this entanglement detection 
scheme in optical lattices and magnetic microtraps fol- 
lows four steps: (i) Creation of two identical copies of the 
entangled state pi23...jv : Each of the two rows of bosons 
shown in Fig. ^ is realized by a ID chain of entangled 
atoms. The entanglement can e.g. be created by spin 
selective movement and controlled interactions between 
atoms as described in |2j or by entangling beam split- 
ters as investigated in [ll| . We assume that any hopping 
of atoms between the lattice sites is initially turned off 
and that the two chains consist of exactly one atom per 
lattice site 12]. (ii) Implementation of the pairwise BS: 
This is achieved by decreasing the potential barrier be- 
tween the two rows of atoms. In an optical lattice one can 
decrease the corresponding laser intensities 0] while in 
an array of magnetic microtraps electric/magnetic fields 
can be switched to change the barrier height [13]. The 
dynamics after lowering the potential barrier is described 
by the Hamiltonian H = Hbs + #int where (c.f. [l^ ) 

N 

H BS = £ -J(a?K$ + b^b^ + h.c.) 
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^bMbPhPb® + Urtb^bPaWaP. (5) 

Here Hbs describes vertical hopping of particles between 
the two rows with hopping matrix element J [l^j and 
Hint gives the on-site interaction of two particles in 
a lattice site with internal state dependent interaction 
strengths U a , Ub and U a b- For simplicity we assume 
that the interaction terms iJi n t can be neglected while 
the hopping is turned on, i.e. J 3> U , and we assume 
U a = Ub = U a b = U 0]. Turning on the Hamiltonian 
Hbs foi the specific time T = tt/(AJ) implements the N 
pairwise BS. (iii) Acquisition of a relative phase between 
the symmetric and antisymmetric parts of the wave func- 
tion: After implementing the BS we let the system evolve 
according to H mt for time r. This introduces a phase 
9 = Ut in each doubly occupied lattice site while it has 
no influence on singly occupied lattice sites. Recently it 
was demonstrated in an interference experiment [ltij that 
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this phase 9 allows the double occupancy sites to be dis- 
tinguished from single occupancy ones, (iv) Turning off 
the lattice and measuring the resulting interference pat- 
tern 0| : After time r the particles are released from the 
trap such that their wave function dominantly spreads 
along the vertical direction x (see Fig. QJ. The density 
profile resulting from the pair of atoms j in state p'j will 



exhibit interference terms dependent on c 



(J) M 



and 6, 



so by varying the interaction phase 9, and noting that 



Pi 



1 — P J , we can determine |c 



Cj)|2 



For N pairs of 



atoms the density profile will depend on |cj^| 2 and |c_ | 2 
as well as on correlations between the density profiles of 
different pairs of atoms, according to Eq. @. Measur- 
ing the density profile for the N pairs of atoms allows 
us to determine the different joint probabilities P±±...±, 
so by solving Eq. (@J we can detect multipartite entan- 
gled states which violate Eq. (0. We note that both the 
creation of the two copies of P123...N and the network 
can be implemented with current experimental technol- 
ogy and do not introduce any novel or unknown sources 
of imperfections. 

We now focus on the predictions of the entanglement 
detection network for multipartite entangled states re- 
cently generated in optical lattices by letting atoms in 
the chain interact with their nearest-neighbors |2|. The 
state generated for an arbitrary number of atoms N is 
given by (l + e^)/2 | 000...0) + (l-e^)/2 | Cluster State), 
where <p is the phase arising from the controlled inter- 
actions during the entanglement creation procedure and 
the precise definition of the cluster state is given in [l7j . 
The entangled state of N atoms generated by the pro- 
cess described in Q will violate the inequalities Eq. 
ViV. For N = 2 tracing out one atom yields a viola- 
tion V of Eq. (J2J for <fi ^ {0,7r, 2tt}, i.e. our network 
detects all entangled states and is thus more powerful 
than Bell-CHSH inequalities which are only violated for 
4> G [if, tt — ip] U [tt + <p, 2n — ip] with tp sa 0.7. When 
starting with N = 3 tracing marginal atoms yields vi- 
olations V as shown by the dashed and solid curves in 
Fig. EI a), respectively. Again, the network detects and 
distinguishes all multipartite entangled states created in 
this case. Also for N > 3 multipartite entanglement is 
detected by tracing the marginal atoms. However, trac- 
ing further atoms does not yield additional violations of 
Eq. (J2J) as indicated by the grey line in Fig.^a), but the 
purity of the further reduced states remains constant. 
This is due to the characteristic of cluster states b eing 
resilient to measurements on individual subsystems |17| . 
In fact, the knowledge of the cluster states' purities can 
be further exploited to identify defects (empty sites) in 
the lattice and errors in the entanglement creation pro- 
cess |Ti| . 

Let us also note that any pure maximally entan- 
gled multipartite state can unequivocally be identified 
through its unique property of being pure while all re- 
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FIG. 2: In Fig. 2(a), we plot the violation V of the in- 
equalities Eq. ©, Vi = Tr(p? 33 ) - Tr(pi 2 ) (dashed), V 2 = 



Tr (pl 2 ) - Tr (pi) (grey) and V 3 



Tr(p? 2 ) -Tr(p|) (solid), 



as a function of the phase <f>, for N = 3 atoms. Whenever 
V > 0, entanglement is detected by our network. In Fig. 2(b), 
we plot the purity IIjv_ m for m = 1 (solid black), m = 7 
(dashed black), m — 14 (solid grey) and m — 20 (dashed 
grey), as a function of e, for N = 300 atoms. 



duced density operators have purity 1/2. In optical lat- 
tices such maximally entangled generalized GHZ states 
that could be tested with the network can e.g. be created 
by beam splitter setups in ID atomic pipelines as shown 
in 

Our network can also be used to study superpositions 
of distinct quantum macroscopic states which are of great 
importance for the better understanding of fundamental 
aspects of quantum theory 0, 0] . There have been sev- 
eral proposals on how to create macroscopic superposi- 
tions in systems ranging from superconductors |2fJ , Bose- 
Einstein condensates (BECs) |21| to opto-mechanical se- 
tups j^. In the case of BECs, the macroscopic super- 
positions are multipartite entangled states of the form 



(16 



+ 



(6) 



and we define a pa- 
2 . Recently, 



where if = 2+ (oh | <p 2 ) -t-\<P2|¥>i 
rameter e by the overlap e 2 = 1 — | ( 4>i \ <p2 
a measure based on e for the effective size S of such su- 
perpositions of distinct macroscopic quantum states was 
introduced jl^ . It compares states of the form | ip) with 
generalized GHZ states of N atoms (| 0)® W + 1 \)® N )/y/2, 
where e = 1 for a generalized GHZ state. The effective 
size S of the state | ip) is given by S = Ne 2 |2^. 

We can determine S from the measurement of the pu- 
rity of any reduced density operator of Eq. © . We derive 
an explicit formula for the purity H]y- m = Tr (p%_ m ), 
where PN-m is the density operator p^ = \ ip) (ip | re- 
duced by m subsystems. We find 



n 



N—m 



1 + 7" 1 + 7 W + 4j N / 2 + 7 
2(1 + 7^/2)2 



N—m 



(7) 



with 7 = 1 - e 2 = I (0x I 2 ) | 2 . 

Suppose we create two identical BECs, each in state 
Pn, wait for a time t c to let their density operators be 
inelastically reduced via single particle loss processes to 
PN-m <8 PN-m', and then let the two BECs go through 
a BS like transformation. As an aside we note that the 
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reduced density operators emerging from multi particle 
collisions not only depend on e, but also on \<f>i) , |0 2 ) 
and thus could be used to gain further insight into the 
properties of the state. The BS can be implemented ei- 
ther through collisional interactions between the atoms 
in two arms of a spatial interferometer |2^. or by first 
turning both BECs into Mott insulator states ^] trap- 
ping them in an optical lattice and then switching on 
Hbs- We only consider the latter method since it cor- 
responds more directly to the situation of Fig. ^ The 
loss processes which reduce the density operators pn are 
stochastic so in general m ^ m! which means that only 
N — n, where n — max{m, m'}, pairs of atoms will un- 
dergo pairwise BS in the lattice. Since only density pro- 
files of pairs of vertical sites with two atoms contribute to 
the interference pattern, measuring the collective density 
profile will determine Tr (p%_ n )- Plots of different Hn-u 
for an initial number of N — 300 atoms as a function of 
e are presented in Fig. Efb). The dependence of these 
curves on N is very weak but for constant e the values of 
IIjV_ n quickly tend towards 1/2 as n increases. Therefore 
from measuring the density profile the determination of 
e from iljv-n is best done for small n ~ 15. For a given 
particle loss rate the average value of n after time t c will 
be known and e can be found by averaging over several 
runs of the experiment performed under identical initial 
conditions. We note that this measurement is consid- 
erably simpler than those in the previous entanglement 
detection schemes, since we do not require the ability to 
distinguish between individual pairs of bosons but only 
need to find the overall probability of projecting on the 
symmetric and antisymmetric subspaces. If the experi- 
mental setup allows to determine the number of pairwise 
beam splitters N — n, the measurement can be performed 
in one run and inelastic processes are not necessary if one 
can measure the collective density profile associated with 
a subset of the pairs of atoms. 

We have presented and investigated a simple quantum 
network that detects multipartite entanglement, requir- 
ing only two identical copies of the quantum state and 
pairwise BS between the constituents of each copy. We 
have shown how the network can be implemented in opti- 
cal lattices and magnetic microtraps, using current tech- 
nology. As examples of its power we have applied the net- 
work to detect entanglement and imperfections in cluster 
states and shown that it also can be used to characterize 
macroscopic superposition states. 
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